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Abstract. We investigate the form factors and decay rates of exclusive ^0 
04 , semileptonic meson decays using the constituent quark model based on the light-front 

quantization. Our model is constrained by the variational principle for the linear plus 
. , Coulomb interaction motivated by QCD. Our numerical results are in a good agreement 

' with the available experimental data. 

^ ■ One of the distinctive advantages in the hght-front approach is the well- 

' estabhshed formulation of various form factor calculations using the well-known 

Drell-Yan-West (^'^=0) frame [1]. In g+=0 frame, only parton-number-conserving 



X 



I Fock state (valence) contribution is needed when the "good" components of the 

current, J"*" and J±={Jx, Jy), are used [2]. For example, only the valence diagram 
shown in Fig. 1(a) is used in the light-front quark model (LFQM) analysis of 
r-| ! spacelike meson form factors. Successful LFQM description of various hadron form 

factors can be found in the literatures [3-5]. 

However, the timelike (g^ > 0) form factor analysis in the LFQM has been 
hindered by the fact that q~^=0 frame is defined only in the spacelike region 
{q^=q^q^ — g^. < 0). While the g'*"7^0 frame can be used in principle to compute 
the timelike form factors, it is inevitable (if g^T^O) to encounter the nonvalence 
diagram arising from the quark-antiquark pair creation (so called "Z-graph"). For 
example, the nonvalence diagram in the case of semileptonic meson decays is shown 
in Fig. 1(b). The main source of the difficulty, however, in calculating the non- 
valence diagram(see Fig. 1(b)) is the lack of information on the black blob which 
should contrast with the white blob representing the usual light-front valence wave 
function. In fact, we noticed [2] that the omission of nonvalence contribution leads 
to a large deviation from the full results. 

In this paper, we circumvent this problem by calculating the semileptonic pro- 
cesses in g+=0 frame and then analytically continuing to the timelike region. The 
g"'"=0 frame is useful because only valence contributions are needed. However, one 
needs to calculate the component of the current other than J"*" to obtain the form 



p. p; p; p; 

(a) (b) 
FIGURE 1. The LFQM description of a electroweak meson form factor: (a) the usual light-front 
valence diagram and (b) the nonvalence (pair-creation) diagram. The vertical dashed line in 

(b) indicates the energy-denominator for the nonvalence contributions. While the white blob 
represents the usual light-front valence wave function, the modeling of black blob has not yet 
been made. 

factor Since J~ is not free from the zero- mode contributions even in q^—Q 

frame [6,7], we use instead of J~ to obtain /_. 

The key idea in our LFQM [5] for mesons is to treat the radial wave function 
as a trial function for the variational principle to the QCD-motivated Hamiltonian 
saturating the Fock state expansion by the constituent quark and antiquark. The 
spin-orbit wave function is uniquely determined by the Melosh transformation. 
We take the QCD-motivated effective Hamiltonian as the well-known linear plus 
Coulomb interaction given by 



H,, ^Ho + = Jml + k^ + J ml + k^ + V, 



qq — -'-'U ^ " qq — y '"-q ^ ro ^ y '"-g -r n -r v qq, 

where 

y„=y. + y,,, = <. + ..-_ + ^vVc... (2) 

We take the Gaussian radial wave function 4>{k'^) = exp(— A;^/2/5^) as our trial 
wave function to minimize the central Hamiltonian [5]. Since the string tension 
6=0.18 GeV^ and the constituent u and d quark masses niu— ma— 0.22 GeV are 
rather well known from other quark model analyses commensurate with Regge 
phenomenology [8], wc take them as our input parameters. The model parameters 
of a, K, and are determined by the variational principle using the masses of p 
and TT [5,9]. It is very important to note that all other model parameters such as 
rric, mb, Puc, Pub, etc. are then uniquely determined by our variational principle as 
shown in [9]. More detailed procedure of determining the model parameters and 
ground state meson mass spectra can be found in [5,9]. 

Our predictions of the ground state meson mass spectra [9] are in a good agree- 
ment with the available experimental data. Furthermore, our model predicts the 
two unmeasured mass spectra of ^So{bb) and 5*1 (6s) systems as Mf,f;=9657 MeV 
and Mbs=5424 MeV, respectively. Our values of the decay constants [9] are also in 



a good agreement with the results of lattice QCD [10] anticipating future accurate 
experimental data. 

The matrix clement of the current j'^ = q2l^Qi for ^~{Qiq) 0~(g2?) decays 
can be parametrized in terms of two hadronic form factors as follows 



(Pi + P2r - ^^^^q^ 



2 



+ Mq') ' , ' q', (3) 



where q^ = (Pi — P2)'^ is the four-momentum transfer to the leptons and mf < 
q^ < (Ml — M2Y. The form factors /+ and /o are related to the exchange of l" 
and 0"*", respectively, and satisfy the following relations: 



M'(-Ml 



/+(0) = /o(0), fo{q') = Uiq') + -TTo^f-il')- (4) 



In the LFQM calculations presented in Ref. [11], the ^''''t^O frame has been used to 
calculate the semileptonic decays in the timelike region. However, when the q~^^0 
frame is used, the inclusion of the nonvalence contributions arising from quark- 
antiquark pair creation (see Fig. 1(b)) is inevitable and this inclusion may be very 
important for light-to-light and heavy-to-light decays. Nevertheless, the previous 
analyses [11] considered only valence contributions in ^'"''7^0 frame neglecting non- 
valence contributions. In this work, we circumvent this problem by calculating the 
processes in q'^=0 frame and analytically continuing to the timelike region. The 
q'^=0 frame is useful because only valence contributions are needed. However, one 
needs to calculate the component of the current other than to obtain the form 
factor f-{q^). Since J~ is not free from the zero- mode contributions even in gf+=0 
frame [6], we use instead of J~ to obtain /_. In the q'^—Q frame, we obtain 
the form factors f+{q'^) and f-{q'^) using the matrix element of the and "_L"- 
components of the current, J'^, respectively, and then analytically continue to the 
timelike q^ > region by changing q± to iq± in the form factors. 

Light-to-light decays: For Ki^ decays, the three form factor parameters, i.e., 
A+, Ao and ^a, have been measured using the following linear parametrization [12]: 

f±{q')-f±{q'-m]){l + X±^), (5) 

where A±,o is the slope of /±,o evaluated at = mf and — f-/ f+\q^=mf- Our 
predictions of the parameters for Ki^ decays in ^+=0 frame, i.e., /+(0), A+, Aq, 
{r'^)K-K = 6/|(0)//+(0) = 6A+/M2+, and ^a = /-//+|g2=^2, are summarized in Ta- 
ble 1. Our result for the form factor /+ at zero momentum transfer, /+(0) = 0.962, 
is consistent with the AdemoUo-Gatto theorem [13] and also in a good agreement 
with the result of chiral perturbation theory [14], /+(0) = 0.961±0.008. Our results 
for other observables such as A+, ^a, and T{Kis) are overall in a good agreement 



TABLE 1. Model predictions for the parameters of Ki^ decay form factors obtained from 
q~^=0 frame. The charge radius r^x is obtained by (r^),ri<' = Qf+{Q^ = 0)//+(0). For 

comparison, we include the results (in square brackets) of the valence contribution in q'^y^O 
frame. The CKM matrix used in the calculation of the decay width (in units of 10^ s~^) is 
\Vus\ = 0.2205 ±0.0018 [12]. 



Observables Our model Other models Experiment 



/+(0) 


0.962[0.962] 


0.961 ± 0.008 [^^1, 0.952 [151,0.93 I^^l 




A+ 


0.026[0.083] 


0.028 [1^1,0.019 [1^1 


0.0286 ± 0.0022[if+3] 
0.0300 ±0.0016[/s:03] 


Ao 


-0.009[-0.017] 


0.0026 [i-'^l, -0.005 [1^1 


0.004 ± 0m7[K+J\ 
0.025 ± 0.006 [/fOg] 




-0.41[-1.10] 


-0.28 [151,-0.28 [15] 


-0.35±0.15[K+3] 








-0.11±0.09[if03] 


(r)„K(fm) 


0.56 


0.57 [15], 0.48 [16] 




7.36 ±0.12 




7.7±0.5ra 



with the experimental data [12]. For comparison, we also include the results (in 
square brackets of the second column of Table 1) of /+(0), A_|_, Aq, and obtained 
from the valence contribution in frame. Even though the form factor /+(0) 

in q^^O frame is free from the nonvalence contributions, its derivative at q'^=0, 
i.e., A+, receives the nonvalence contributions. Moreover, the form factor f-{q^) 
in ^'''t^O frame is not immune to the nonvalence contributions even at ^^=0 [6]. 
Unless one includes the nonvalence contributions in the q^^O frame, one cannot 
really obtain reliable predictions for the observables such as A_|_, Aq and for Kis 
decays. 

In Fig. 2, we show the form factors /+ obtained from both q'^—0 and q'^^0 
frames for < < (M^- - M^)^ region. As one can see in Fig. 2, the form factor 
f+{q^) obtained from g+=0 frame (solid lines) [2] appears to be linear functions of 
q^ justifying Eq. (5) usually employed in the analysis of experimental data [12]. 
Note, however, that the f+{q^) obtained from only valence contribution in q'^^Q 
frame (dotted lines) does not exhibit the same behavior. 

Heavy -to -light (heavy) decays: Our predicted decay rates ioi D ^ K and 
D ^ TT are V{D^ K-f+u^) = 8.36|V;,p x lO^^ pg-i and T{D^ vr-e+z/e) = 
0.113|VcdP ps"S respectively. Using \Vcs\ = 1.04 ± 0.16 and = 0.224 ± 0.016 
[12], we obtain the branching ratio of Br{D ^ K) = (3.75 ± 1.16)% and Br{D 
tt) = (2.36±0.34) X lO'^, while the experimental data are (3.66±0.18)% for D ^ X 
and (3.9t?:? ± 0.4) x lO'^. Also, our predicted decay rates for B ^ tt and B ^ D 
are T{B° n-i+ue) = 8.16|K6p ps"\ r(fi° ^ D-f+Ue) = 9.39\Vcb\'^, respectively. 
Using \Vub\ = (3.3 ± 0.4 ± 0.7) x lO^^ and \Vbc\ = 0.0395 ± 0.003 [12], we obtain 
Br{B ^ tt) = (1.40 ±0.34) x 10"^ and Br{B ^ D) ^ (2.28 ± 0.20)%. Our results 
are quite comparable with the recent experimental data [12], (1.8 ± 0.6) x 10~^ for 
B ^ 71 and (2.00 ± 0.25)% for B ^ D, within the given error range. 

In the heavy quark limit Mi(2)^oo, the form factor f+{q^) is reduced to the 
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FIGURE 2. The form factors f+{<f') for the K w transition in timehke momentum transfer 
> 0. The soUd and dotted lines are the results from the g+=0 and q~^^0 frames[2], respec- 
tively. The differences of the results between the two frames are the measure of the nonvalence 
contributions from g'+^O frame. 

universal Isgur-Wise (IW) function, C{vi-V2) = [2y/Mi Ms/ (Mi + M2) ] /+ (?^) , where 
■^1(2) = -Pi(2)/^i(2)- Our prediction of the slope p^=0.8 of the IW function at the 
zero-recoil point defined as ^{vi-V2) = 1 — p^(f 1 - ^2 — 1) is quite comparable with the 
current world average pavg.=0.66±0.19 [12] extracted from exclusive semileptonic 
B^DiP decay. 

In conclusion, we analyzed the exclusive 0~^0~ semileptonic meson decays using 
the LFQM constrained by the variational principle for the QCD-motivated effective 
Hamiltonian with the well-known linear plus Coulomb interaction [5,9]. The form 
factors f± are obtained in q'^=0 frame and then analytically continued to the 
timelike region by changing q± to iq± in the form factors. The matrix element 
of the "_L" component of the current is used to obtain the form factor /_. Our 
model provided overall a good agreement with the available experimental data and 
the lattice QCD results for the transition form factors and branching ratios of the 
0~ 0~ semileptonic meson decays. Also it rendered a large number of predictions 
to the heavy meson mass spectra and decay constants [9] . We think that the success 
of our model hinges on the advantage of light-front quantization realized by the 
rational energy-momentum dispersion relation. It is crucial to calculate the "good" 
components of the current in the reference frame which deletes the complication 
from the nonvalence Z-graph contribution. The present work broadens the utility 
of the standard hght-front frame a la Drell-Yan-West to the timelike form factor 
calculation. We anticipate further stringent tests of our model with more accurate 
data from future experiments and lattice QCD calculations. 
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